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Abstract. We propose a calibrated two-dimensional cellular automaton 
model to simulate pedestrian motion behavior. It is a v max = 4 (3) model 
with exclusion statistics and random shuffled dynamics. The underlying 
regular grid structure results in a direction-dependent behavior, which 
has in particular not been considered within previous approaches. We 
efficiently compensate these grid-caused deficiencies on model level. 

Key words: multi-agent, transition matrix, stochastic approach, cali- 
brated model 



1 Introduction 

The different model approaches for microscopic person dynamics are based on the 
particular discipline analogies, ranging from hydro-dynamic models to artificial 
intelligence and multi-agent systems [1]. The complex dynamic human behav- 
ior is induced by individual decisions, which are classified to be of short-range 
(operational) and long-range type (strategic/tactical). The self-organization of 
persons is a further essential characteristic of human behavior. In contrast to 
the social force model [2-4] or the discrete choice model [5, 6] the developed mo- 
tion model [7-9] is based on a stochastic approach to handle the unpredictable 
behavior by individual path deviations. The stochastic motion model is an appro- 
priate and fast method for analysis the dynamic pedestrian behavior. However, 
to derive valid results several simulation runs (>100) have to be performed. The 
focus concentrates on the evaluation of application oriented simulation scenar- 
ios instead of the characteristics of individual interactions or specific pedestrian 
trajectories. 

2 Stochastic Motion Model 

The presented motion model is based on a stochastic approach [10], which is 
comparable to a common cellular automaton. It utilizes a regular grid structure. 
In contrast to the cellular automaton, the new model is developed on the basis 
of a fundamental paradigm shift: instead of changing the cell status depending 
on the status of its surrounding cells (neighbors), the agent is able to move 



over the regular lattice and to enter those cells, which are not occupied by 
other agents or obstacles (e.g. walls). To describe the motion behavior of an 
agent, the motion vector is separated into a desired motion direction and a 
transversal deviation [10]. Using the spatially discrete grid structure and defining 
three transition states (forward | stop | backward or left | on track | right) the 
normalized transition probability (p) into these states is generally defined by the 
following equations. 

p + = — (a 2 + fj 2 + /i) , for forward or left 
p° = 1 — (cr 2 + /t 2 ) , for stop or on track (1) 
p~ = — (cr 2 + fi 2 — fx) , for backward or right 

In the case of the desired motion direction, /j, denotes the desired speed 
and cr 2 the corresponding variance. If the transversal deviation is concerned, 
fi is the average and a 2 is the range of the fluctuations. Considering a sym- 
metric transversal deviation (//deviation = 0) an d a connection of desired speed 
and the corresponding variance (no step backward p~ = 0, so that c 2 pocd = 
Aispecd(l — Mspeed), the above equations are simplified to the following equations 
for the desired motion direction (2) and for the transversal motion direction (3). 



p fOTWard = Mspeed I p St ° P = 1 ~ Mspeed (2) 

left, right _ ± 2 | on track _ -i _ 2 (o\ 

P — 2 deviation I P — 1 "deviation \° ) 

Finally, the motion components are combined to a 3x3 transition matrix 
(My) as shown in the following fig. 1. The emphasized cell (marked gray at the 
figure) contains the transition probability of moving forward without transversal 
deviations. In fact, the transition matrix possesses a two-dimensional character- 
istic, but it only defines an one-dimensional transition considering a transversal 
deviation (1.5-dimensional). 
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Fig. 1. Generation of the transition matrix due to combination of desired motion 
speed and transversal motion deviation. 



To allow for a three-dimensional agent motion behavior, two independent mo- 
tion directions are needed. Based on the developed horizontal transition matrix 
(M) a diagonal transition matrix (M) is derived by re-indexing the horizontal 
matrix (turning the matrix by ^, [10]). The motion direction (a) is integrated 
into the stochastic model by superpose these matrices with A. 

!(1-A)M + AM, A = tana, < a < f 

(4) 
-±(1-X)M + V2XM , A = tan(f -a),f <a< f 

The rotation of M (4-fold symmetry) allows for determining the entire spec- 
trum of the motion direction. The underlying regular grid structure results in 
a direction dependent behavior (e.g. entering diagonal cells implies walking a 
longer way in comparison to horizontally located cells). Therefore the first model 
modification is to adjust A within the parameters of ? < a < ^ . 

It's obvious, that the stochastic motion model allows for horizontal and di- 
agonal movements (fig. 2. Moore- Neighborhood) . 




Fig. 2. The distance metric d results from the Moore-Neighborhood. 

The shortest distance between two points (start, destination) is given by the 
metric d (5) which differs from the defined p-norms (Manhattan, Euclidean and 
Chebychcv norm). 

d = | An - Ay\ + v / 2min(Ax, Ay) (5) 



3 Model Constraints, Improvements and Validation 

Due to the utilization of a regular grid structure, the transition matrix does 
not fulfill the criteria of independent agent motion behavior. So, the speed and 
the variance of the agent depend on the agent motion direction. If the agent 
enters diagonal cells his walking distance is longer (approx. 41%) in comparison 
to the use of horizontally located cells. This model constraint is equivalent to a 



significant higher motion speed depending on the direction of motion. Algorithms 
to compensate this grid-based speed effect can be found at [7, 11]. A detailed 
model analysis points out that the expected value of the transition matrix /xm 
(6), defined by cell based transition probability My and the relative location 
e,j, differs from /x S pccd, which is specified in the motion model. 

W' = EE e y M ^ ' with eij= (j) ( 6 ) 

i j ^ ' 

Fig. 3 points out the correlation of speed and motion angle. With increasing a 
the corresponding speed (/J. S peed) increases as well, whereas the stepwise change of 
the transversal deviation (cj eviation ) mitigates the a-depending characteristics. 




Fig. 3. Model deficiency due to motion angle a depending characteristic of speed. 



The stochastic motion model allows for absolute cell transitions (My = 1 
with CTj eviation = 0) only at motion angle of a = and a = f . If the agents 
choose another angle he always has to choose between two cells at least. Because 
of the superposition of the horizontal (M) and diagonal (M) transition matrices, 
a model immanent variance from desired agent motion direction occurs, even if 
the model defines "deviation = 0). The characteristic of this overall motion model 
variance is shown in fig. 4. 

The angle depended variance of motion implies to major issues. The model 
shows an immanent motion deviation, which is not considered in previous equa- 
tions and the different variances lead to different avoiding behavior. Using the 
parameter set ^t spcc d = 1 and ^deviation = as an example, at a = no variance 
is allowed by the model and an agent cannot move if the chosen cell is blocked. 
Using the same scenario with a different motion angle (e.g. a = ^), the agent gets 
the probability of approx. 20 % to pass the blocked cell. To ensure homogeneous 
variance, an appropriate compensation on model level is needed. The expected 




Fig. 4. Model deficiency due to motion angle depending characteristic of vari- 
ance. 



value fiM of the matrix depends on /J, sp eed and cr^ oviation whereas the parameters 
are directly coupled. For each parameter set a specific characteristic over the 
angle has to be calculated. The following fig. 5 shows these characteristics. 

Further model investigations point out that the model has to be extended to 
reproduce the representative shape of the fundamental diagram [12, 13]. Using 
the transition matrix, an agent is able to react to the status of the adjacent cells 
(empty, occupied). The fundamental diagram indicates an interaction range of 
about 1.3 m, because the speed of an agent starts to decrease if the density 
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Fig. 5. Functional characteristics of the parameter compensation for the param- 
eter set /i spccd = 0.9 and cr^ eviation = 0.2. 



relations /o/p max reaches a level of 10 % (considering an agent with a dimension 
of 0.4x0.4 m and a maximum density of p max = 6.25 Person/m 2 ). If the agent 
moves three/four steps at once, he will be able to interact with distant agent 
and the developed model is found to reproduce the characteristic shape of the 
fundamental diagram (fig. 6). Therefore the motion model has to provide the 
following agent properties: 

— Always move, no waiting (occupied cells increase transition probability of 
the other matrix cells). 

— Move four steps at once and decrease the steps depending on agent density, 
at least in the case of p/p ma x > 0.6 the number of steps should be reduced 
from four to three to fit the shape. 

— Agent leaves a trace, at each time step all entered cells will temporarily 
blocked. 
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Fig. 6. Fundamental diagram of stochastic motion model using different agent 
operating distances (3-4 steps at once) depending on agent density. 

The stochastic model meets all criteria for a scientifically reliable motion 
model. It exhibits the absence of significant model-caused limitations and re- 
produces all common self-organizing effects (e.g. row formation or oscillation). 
Besides the operational motion definition by the stochastic transition matrix, 
strategic/tactical motion components are taken into account as well. The stochas- 
tic model allows for the reaction of the agent to objects/agents at immediate 
vicinity. It additional provides the capability of considering distant constellation 
of agents (jam) and potentially blocked bottlenecks. 

4 Summary and Outlook 

Model specific parameter corrections ensure that the motion vector is equal 
to the expected value of the corresponding transition matrix. This issue has 



in particular not been considered within previous approaches. The calibrated 
motion model is thus the hrst approach, which allows for a specific stochastic 
description of agent movements without model restrictions. 

The passenger related evaluations and simulations of dispatch processes at 
Dresden Airport exemplarily show that the developed stochastic motion model 
is able to reproduce the behavior of passengers in an appropriate way [8, 9]. 
Therefore the developed motion model is implemented in a corresponding ap- 
plication environment which allows for a various application, e.g. investigations 
regarding to group dynamic behavior or route planning in the airport terminal 
focused on normal operations and emergency cases. 



Bibliography 



[1] A. Schadschneider, W. Klingsch, H. Kluepfel, T. Krctz, C. Rogsch, and 

A. Seyfried, Evacuation Dynamics: Empirical Results, Modeling and 

Applications, ch. Evacuation Dynamics: Empirical Results, Modeling and 

Applications, pp. 3142-3176. Springer, 2009. 
[2] D. Helbing and P. Molnar, "Social force model for pedestrian 

dynamics", Phys. Rev. E 51 (May, 1995) 4282-4286. 
[3] A. Johansson, D. Helbing, and P. Shukla, "Specification of the social force 

pedestrian model by evolutionary adjustment to video tracking data" , 

Advances in Complex Systems 10 (2007) no. 4, 271-288. 
[4] M. Moussaid, N. Perozo, S. Gamier, D. Helbing, and G. Theraulaz, "The 

walking behaviour of pedestrian social groups and its impact on crowd 

dynamics", PLoS One 5 (2010) no. 4, . 
[5] M. Bicrlaire, G. Antonini, and M. Weber, "Behavioral dynamics for 

pedestrians", in Proceedings of the 10th International Conference on 

Travel Behavior Research. Lucerne, 2003. 
[6] T. Robin, G. Antonini, M. Bierlaire, and J. Cruz, "Specification, 

estimation and validation of a pedestrian walking behavior model" , 

Transportation Research Part B: Methodological 43 (2009) no. 1, 36-56. 
[7] M. Schultz, S. Lehmann, and H. Fricke, Pedestrian and Evacuation 

Dynamics 2005, ch. A discrete microscopic model for pedestrian dynamics 

to manage emergency situations in airport terminals, pp. 369 - 375. 

Springer, Berlin, 2007. 
[8] M. Schultz, Entwicklung eines individuenbasierten Modells zur Abbildung 

des Bewegungsverhaltens von Passagieren im Flughafenterminal. PhD 

thesis, Technische Universitat Dresden, 2010. 
[9] M. Schultz, C. Schulz, and H. Fricke, Pedestrian and Evacuation 

Dynamics 2008, ch. Passenger Dynamics at Airport Terminal 

Environment, pp. 381 - 396. Springer, Berlin, 2010. 
[10] C. Burstedde, K. Klauck, A. Schadschneider, and J. Zittartz, "Simulation 

of pedestrian dynamics using a 2-dimensional cellular automaton" , 

Physica A 295 (2001) 507-525. 
[11] T. Kretz and M. Schrcckcnbcrg, Pedestrian and Evacuation Dynamics 

2005, ch. Moore and more and symmetry, pp. 297 - 308. Springer, Berlin, 

2007. 

[12] U. Wcidmann, "Transporttcchnik der Fufiganger", Schriftenreihe des 

Institut fur Verkehrsplanung, Transporttechnik, Strassen- und 

Eisenbahnbau 90 (1992) . 
[13] A. Seyfried, B. Steffen, W. Klingsch, and M. Boltes, "The fundamental 

diagram of pedestrian movement revisited", J.STAT.MECH. (2005) 

P10002. 



